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HIGH ENERGY NEUTRINO AND GAMMA-RAY TRANSIENTS 
FROM RELATIVISTIC SUPERNOVA SHOCK BREAKOUTS 
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ABSTRACT 

Relativistic shocks that accompany supernovae (SNe) produce X-ray burst emissions as they break 
out in the dense circumstehar medium around the progenitors. This phenomenon is sometimes as- 
sociated with pecuhar low- luminosity gamma-ray bursts (LL GRBs). Here, we investigate the high 
energy neutrino and gamma-ray counterparts of such a class of SNe. Just beyond the shock breakout 
radius, particle acceleration in the collisionless shock starts to operate in the presence of breakout 
photons. We show that protons may be accelerated to sufficiently high energies and produce high 
energy neutrinos and gamma rays via the photomeson interaction. These neutrinos and gamma rays 
may be detectable from < 10 Mpc away by IceCube/KM3Net as multi-TcV transients almost simul- 
taneously with the X-ray burst emission, and even from < 100 Mpc away with follow-up observations 
by CTA using a wide-field sky monitor like Swift as a trigger. A statistical technique using a stacking 
approach could also be possible for the detection, with the aid of the SN optical/infrared counterparts. 
Such multi-messenger observations offer the possibility to probe the transition of relativistic shocks 
from radiation-mediated to collisionless ones, and would also constrain the mechanisms of particle 
acceleration and emission in LL GRBs. 

Subject headings: acceleration of particles — neutrinos — shock waves — supernovae: general 



[ 1. INTRODUCTION 

There is increasing evidence for a class of super- 
novae (SNe) which are accompanied by relativistic shocks 
(^sh ^ 0.1). These are initially radiation mediated while 
going through optically t hick gas, i.e., the stellar or 
a circumstehar envelope (jWeaveii Il976[) . The typical 
photon energy could be as high as 10-100 keV within 
a few Thomson mean free paths nea r the shock (e.g., 
iKatz et al.l 120101 : iBudnik et aLll2010[ ). Once the opti- 
cal depth becomes small enough that the photon dif- 
fusion velocity is similar to the shock velocity, photons 
begin to escape and are observed as the shock break- 
out emission. The energy spectrum is expected to be 
a combination of non-thermal and quas i- thermal (e.g., 
Wang et al.l 120071: IWaxman et al.l l2007t iNakar fc Sad 
2010i: iChevalier fc IrwinI l2012f ). Such X-ray breakouts 



have now been confirm ed in some cases ([Soderberg et al.l 
[200l[OfekeraII[20T2h . 

An interesting application of relativistic shock break- 
outs is to low luminosity gamma-ray bursts (LL GRBs) 
associated with SNe, s uch as GRB0602 1 8/SN 2006a.i 
(|Campana et al.l 120061: ISoderberg et al.l I2006D and 
GRB100316D/SN20 10bh (jChornock etall [20T0I: 
IStarling et al.l 1201 ID . These GRBs share peculiar 
characteristics compared to typical long GRBs: smaller 
luminosities (Liso,7 ~ lO"'^ erg/s), softer peak energies 
(cpeak ^ 1-10 keV), and longer durations {t^ ~ 1000 sec). 
The relativistic shock breakout model can explain these 
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characteristics (jSoderberg et al.l I2006D although a 
more det ailed modeling is ne cessary for the energy 
spectrum (jGhisellini et al.ll2007[) . The host galaxi es are 
both nearby (z = 0.0331 for GRB060218 (Mirabal et al 
20061) and z = 0.059 for GRB100316D ^Vergani et al 



2OIOD ). The dimness and the proximity indicate that 



there might be a considerable number of sub-threshold 
events resulting in a relatively high local rate of some 
Rj,t(z = 0) < 500 G pc~^vr-i (IGuetta fc Delia VaUd 
120071: iLiang et al.ll2007t) . Hence LL GRBs may be one of 
the most promising targets for ongo ing and upcomin g 
multi- messenger astronomy searches (jAndo et al.ll2012f ). 
e.g. with specially designed facilities such as AMOP>0. 

In general, shocks would become collisionless beyond 
the break out radii, allowing charged particles to be ac- 
celerated (|Murase et al.ll2011l : lKatz et al.ll2011| ). The ac- 
celerated protons can then produce mesons by the pho- 
tomeson interaction or the inelastic pp reaction, result- 
ing in neutrinos and gamma rays as decay products of 
the mesons. This situation is analogous to the photo- 
spheric scenario of GRB prornpt em ission (iMuras 3 120081: 
IWang fc Dal 120091 : iGao et al.ll2012[ ). By detecting these 
signals simultaneously with the breakout emission, one 
can probe the evolution of shocks converting from the 
radiation mediated to the collisionless regime, and inves- 
tigate the baryon acceleration there. 

In this Letter, we evaluate the high energy neutrino 
and gamma-ray emission from relativistic shocks break- 
ing out. We are primarily concerned with the application 
to LL GRBs associated SNe, but our prescriptions can be 
generally applied to other situations, e.g., ordinary (non- 
GRB) Ibc SN shocks breaking out of the progenitor, by 
appropriately scaling the velocity and the breakout ra- 



^ IMazzali et al.M2006l) and lToma et al.1 I I2007I) proposed strongly 
magnetized neutron stars, instead of black holes, as the ce ntral 
engine to explain the characteristics (see also lFan et al.ll2011l ). 

^ http://amon.gravity.psu.edu/ 
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dius of the shock. 

Our work differs from pr evious studies (jMurase et alj 
l2006t iGupta fc Zhanej[2007D where neutrino counterparts 
of LL GRBs were evaluated based on the relativistic jet 
model (Fi > 1] whic h is also viable (|Toma et all 120071: 
iGhisellini et al.ll2007| ). Although the main microphys- 
ical process, i.e. the photomeson interaction between 
accelerated protons and the observed X-ray photons, 
is the same, the predictions are not the same, due to 
diffe rences in the em ission radius and the Lorentz fac- 
tor (jKatz et al.ll201l[ ). Also, our prediction on the de- 
tectability of the gamma-ray counterpart is more posi- 
tive than the previo us estimate. The pre sent results also 
differ from those of iMurase et al.l (|2011[ ). which consid- 
ered the neutrino and gamma-ray emission from non- 
relativistic shocks (/3sh ^ 0.01) colliding with the dense 
circumstellar material 

2. MODEL SET-UP 

Our main parameters are the shock velocity /3sh and 
the shock breakout radius rsb, where the optical depth 
ahead of the shock becomes as small as « l//3sh- The 
duration of the breakout X-ray emission is comparable 
to the dynamical time scale. 



tdy 



2.0 X 10^ r,b,i3.8/3sh"' sec. (1) 



The intrinsic energy budget can be estimated as 
47r(l/2)prsb'/T/3shV, or 



5.1 X 10"*^ r,b,i3.8^^sh erg 



(2) 



Here Zx and ctt are the mean free path and the cross sec- 
tion of the Thomson scattering, and we use the relation 
at the shock breakout radius, 1^ ~ ~ 'mp/o'TpPsh- 
The mean photon luminosity is then iiso.7 ^ ^i^iao/t'y, 
or 

Liso,7 ^ 2.6 X lO''^ e^rsbaa.s^^sh erg/sec, (3) 

where is the fraction of the energy delivered to ra- 
diation. Wc assume that the prompt X-ray emission 
of GRB060218- and GRB100316D-like events is due to 
the shock breakout emission, and set /?sh = 1 and 
Tsh = 6.0 X 10^'^ cm as the fiducial values which roughly 
reproduce the observed characteristics. 

It is difficult to evaluate the time-dependent energy 
spectrum of the photon field in-situ around the shock 
break out radius (jNakar fc Sari|[20Tot IChevalier fc Irwinl 
[20TI . Here, we substitute instead the observed spec- 
trum which, over most of the prompt emission episode, 
can be fitted approximately by a broken power law, 
dn/de = nb(£/eb)~" for e,„in < e < Eb, and dn/de = 
nh{e/£h)~^ for £b < e < £max- Q We set a = 1.42, 
/3 = 2.48, £b = 16 keV and Lfgo,7 that is the lumi- 
nosity at £b, foh owing the observation of GRB100316D 
()Fan et al.l2011|) . Note that Lf^^ ^ can be related to iiso,7 

via U.y = J'^""''^ e {dn/de) de = iiso,7/47rrsb^c. The low 
energy cut-off is assumed to be Emin = 1-0 keV where 
the observed SED declines sharply. In fact. Swift UVOT 
identified a prompt UV-optical emissio n with a photon 
index harder than 3/2 for GRB060218 (jGhisellini et al.l 

^ We have confirmed that the following arguments and results do 
not change dramatically if we substitute a blackbody like spectrum. 



I2007D. and gave onl y an upper limit for GRB100316D 
(jStarling et al.|[20Tll) . These soft photons can affect the 
opacity of high energy gamma rays, as we show below. 
The high energy cut-off can be set as £max = 1.0 MeV 
since the pair absorption would be crucial above this en- 
ergy. 

The mass density of the circumstellar medium at the 
shock breakout radius can be evaluated as 

4.2 X 10-1^ r,b,i3.8"'/3sh"' g/cml (4) 



crT?'sb/3sh 



If the matter is provided by a stellar wind where p « 
M / Anvvirsh'^ , the mass ejection rate may be as large 
as M ^ 0.03 Ww.9''sb,i3.8/3sh~^ Mq/jv for a wind ve- 



locity 



10^ Ww,9 cm/s. Such massive winds in 



the pre-coUapse phase have been inferred es pecially for 
type Iln SN progenitors (ISmith et al.l l2010|:1Ston et al.l 



[2011 [ChalTdra et al.ll2012allbl : IRoming et al. 1 120121) ■ and 

in some cases, the shock break out X-ray photons have 
been observed ()Ofek et al.ll2012D . 

3. PROTON ACCELERATION 

Beyond the breakout radius, the shock can no longer 
be mediated by the radiation. However, electromagnetic 
fields amplified by plasma instabilities will determine 
the shock structure, i.e. the shock be comes collision- 
less (jMurase^EaDIlSil Eatz^TaDIlQ^ The magnetic 
field strength can be estimated as i? « (87r^BC^7)^^^ ^ 

9.6x104 (eB/0.1)i/V,b43.8-'/'/3sh G, where Cb = Ub/U^ 
represents the magnetic field amplification efficiency. 
A fraction of the protons can be accelerated into a 
non-thermal distribution by the first order Fermi pro- 
cess as they repeatedly cross the shock. The maxi- 
mum energy of the protons is determined by the con- 
dition, iacc < ip, where tacc = rjEp/eBc is the acceler- 
ation time (in the Bohm diffusion limit, rj ^ 27r), and 

tp^ = ^ad^ + tsyn~ + tlC~ + ^BH^ + tpj^ + tpp~ 

is the total energy loss time scale, including the effect of 
the adiabatic expansion (ad), the synchrotron radiation 
(syn), inverse Compton radiation (IC), the Bethe-Heitler 
type interactions (BH), the photo meson interactions 
(p7), and the inelastic pp collisions (pp). Fig IT] shows 
the proton energy versus the various t ime scales, ob - 
tained with the calculation codes used in lMuras^ (j2008| ). 
One can sec that the synchrotron cooling time fgyn = 
(3mp'*c'^/4(TT?7iG^)C^B"^£'p~^ is the most important for 
the proton maximum energy. The condition t^cc ^ ^syn 
appropriately gives 

Ep < 6.3 (ll) (|^)~'^Vsb.i3.8^/^/3.h-^/^ EeV, 

(5) 

and protons can be accelerated up to EcV energies. 
We assume a power law distribution of the protons, 
dN/dEp oc Ep^. The normalization is determined by 
introducing the acceleration efficiency, ecR = Scn/^iso- 
Here Sqr = /J''''"""' Ep{dN/dEp)dEp is the total energy 
of the non-thermal protons. We set Ep mm = IGeV, and 
E„ ,max is given by Eg. ([5]) . 

4. NEUTRINO AND GAMMA-RAY EMISSION 
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Fig. 1. — Acceleration time scale and various cooling time 
scales of protons for a relativistic shock breakout with /3ah = 1-0, 
''sb = 6 X lO^'^cm. For the breakout photons, we assume a broken 
power low with a bolometric luminosity L^^^ ^ = Iff*® erg/sec as 
explained in the text. We also sot = 0.1, and r) = 2-k. 

We consider the neutrinos and the gamma rays from 
the decay of mesons generated by both the photomeson 
production and inelastic pp reaction. In the analytical 
estimate below, we only discuss pions which turn out to 
give a dominant contribution. But the co ntributi o n from 
kaon decay is numerically included as in IMurasd (I200I . 

We can estimate the fraction of energy transferred from 
the non-thermal protons to the pions by the photomeson 
interactions as min[l,/p.^] where f n-y = tnyn/^p.^-. Us- 
ing t he rectangular approximation (jWaxman fc Bahcalll 
|1997[) for a photon spectrum approximated as a broken 
power law, we have 

""{{E^/E^.^r-^ {E^,^<Ep), 

where Ep^^ = 0.5 eeh^^mpC^ ^ 8.8 £b,i6kcV^^ TeV with 
£ ~ 0.34 GeV. The multi-pion production becomes dom- 
inant above « 0. 5 ££mi n~^™pC^ ~ 140 £min kev~^ TcV 
(cf. IMurase et all [2OO80 • We can conclude that a sig- 
nificant fraction of non-thermal protons with energies 
10 TeV ^ Ep < EeV will be converted into pions. 

The inelastic pp cooling time is tpp~^ ~ 
(p/mp)Kpp(TppC. The fraction of energy an incident 
proton loses, /pp = tdyn/^pp, can be evaluated as 

/pp ^ 0.03 /?,h-', (7) 

where we use approximately constant values for the in- 
elasticity Kpp ~ 0.5 and for the cross section cTpp ~ 
4 X 10^^^ cm^, appropriate at high energies. Eq.© indi- 
cates that the inelastic pp collisions can also contribute 
moderately to the pion production as in the case of GRB 
photospheric emissions. 

Neutrino Emission.- Neutrinos are mainly pro- 
duced as decay products of charged pions. One 
can find that the charged pions with E^^ > 

2.7 (CB/0.1)-i/2e^-i/2r,b,i3.8'/'/3sh-' PeV wih lose their 
energy before decaying due to the synchrotron cool- 
ing. Given that the resultant neutrinos have typically 
^ 1/4 of the parent pion energy, one expects TeV- 
PeV neutrinos. The peak fluence from a single su- 
pernova/burst event can be analytically estimated as 



E,^(f. ~ (1/4^^E) X {l/i)mm[lJp^]{Ep^dN/dEp), or 
E.'^'. - (yoS^) /P7^sb,i3.8'/3sh erg/cm^ 

where Z?l is the luminosity distance to the source. 

Fig 12] shows the energy fluences of neu trinos obtaine d 
numerically using the calculation codes of IMurasd (|2008|1 . 
for the same parameters as FiglT] The dashed and 
dotted lines show the contribution from the photome- 
son and inelastic pp interactions, respectively. We 
have verified that contributions from the kaon de- 
cay become important only above ~ 10 PeV. The 
signal is well above the zenith-angle-averaged atmo- 
spheric neutrino background (ANB, dotted-dash lines; 
thick one for ^ 2 x 10'^ sec, and thin one for 
1 day). The number of muon events due to the 
muon neutrinos above TeV energies can be estimated as 
~ 0.2 (ecR/0 2)(AVlOMpc)-Vsh.i.Ts'/3sh using Ice- 
Cube /KM3 net (|Karle fc for the IceCube Collaboration! 
\wWt . \KaMmi^ . Based on our fiducial parameters, Ice- 
Cube/KM3net can marginally detect a nearby source at 

< 10 Mpc, although such events occur very rarely i.e., 

< 0.002 yr-i for a l ocal LL GRB event rate R Lhjz = 
0) ~ 500 Gpc^Vr"^ (IGuetta fc Delia Valid [2007h . 

From Fig 121 one can see that the typical neutrino en- 
ergy in the relativistic shock breakout model is TeV-PeV. 
By comparison, the relativistic jet models of LL GRBs 
predict relatively higher energy PeV-E eV neutrinos 
ijMurase et a"ni200l iGupta fc Zhand 120071 ). This differ- 
ence is mainly because the shock breakout model involves 
a lower Lorentz factor and a stronger cooling of mesons. 
In a relativistic jet, the peak photon energy in the co- 
moving frame is e'^ = £b/rj, where Fj is the Lorentz 
factor of the jet. The typical energy of protons inter- 
acting with photons via the photomeson production is 
Ep' ~ 0.5 ££b'~^'TipC^. The resultant neutrino energy 
will be En ^ 0.05 x Ep'Tj in the observer frame, which 
is 100 (Fj/lO)^ times larger than our model. Thus, high- 
energy neutrino observations can provide clues to the 
emission model of LL GRBs. 

In principle, the shock velocity could be independently 
constrained through the neutrino spectroscopy. From 
Eq. ^ and ^ , both fp~^ and /pp are present irrespective 
of Tsb, and only depend on (3sh- The relative importance 
of photomeson to inelastic pp collisions directly affects 
the neutrino energy spectrum. In the case of relativistic 
shocks with /3sh ^0.1, the spectrum will have a bumpy 
structure like Figl21 On the other hand, slower shocks 
will produce relatively flat spect ra because o f efficient 
inelastic pp interactions (see e.g.. IMurase et ah (2011)). 

Gamma-ray counterparts.- Gamma rays are mainly 
injected by neutral meson decays. Since the neutral 
mesons do not suffer synchrotron cooling, the maxi- 
mum energy of gamma rays can be as high as ~ 10% 
of the parent protons, that is ^ 100 PeV in our fidu- 
cial case. At high energies above ~ MeV, the electron- 
positron pair production can attenuate the gamma-ray 
flux. In the emission region, we can roughly take into 
account the attenuation by w 1/(1 + ^77), where r77 
is th e electron-pos itron pair production optical depth 
fc.g.. lBarinEjl2006l ) . The observed photon spectrum of 
GRB100316D is employed to calculate the optical depth 
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Fig. 2. — Energy fiuences of neutrinos from a relativistic shock 
breakout using the same parameters as Fig[T] We set ecu = 0.2 
and Dl = 10 Mpc. Lines represent a contribution from the pho- 
tomeson production (dashed), the inelastic pp reaction (dotted), 
and the total (solid). The dotted-dashed lines show the zenith- 
angle averaged atmospheric neutrino background (ANB) within a 
circle of dcg for At = 2.0 X 10^ sec (thick) and At = 1 day (thin). 

of the emission region numerically, with a Rayleigh- Jeans 
tail below Emin = 1 kcV. This would be reasonable, since 
the result is not affected much as long as the photon in- 
dex there is harder than 1. Note that the gamma-ray 
attenuation on matter due to the Bethe-Heitler pair pro- 
duction is not important in the case of relativistic shock 
breakouts. The optical depth can be roughly estimated 

as TBH « (l/137)crT(/o/™p)?-sb ~ 0.0074 This can 

be important in the case of non-rela tivistic shocks with 
Psh < 0.01 (e.g.. iMurase et al.|[20lll ). We also take into 
account the attenuation by the e xtragalactic background 
light (EBL: iKlIeiske et al.|[20M ). 

Fig 13] shows the numerically calculated energy spec- 
trum of gamma rays. The thick solid line represents the 
expected flux from a single LL GRB event at 10 Mpc 
and 100 Mpc. The emission duration is set to that of 
the X rays, ~ 2.0 x 10'^ ''sb.is.sAh"^ sec. As a ref- 
erence, we also show the injected spectrum without at- 
tenuation (dashed line) and only including the attenua- 
tion within the emission region (thin solid line) for the 
10 Mpc case. It can be seen that the attenuation of GeV 
^ Ej < 100 TeV gamma rays is mainly due to the pho- 
ton field in the emission region below/around e >1 keV. 
In our case, the attenuation rate decreases with the en- 
ergy because of the Klein-Nishina suppression. On the 
other hand, gamma rays above ^ 100 TeV are mostly 
attenuated by the EEL. In FiglSl we also show the dif- 
ferential sensitivity of CTA for a 5a detection with an ex- 
posure ti me comparable to t^, 0.5hr = 1.8 x 10"^ sec (dot- 
ted line; [Actisiran Honl) . One can see that CTA can 
detect the multi-TeV gamma rays even from 100 Mpc, 
within which the all-sky event rate would be 2 yr~^ 

for i?LL(z = 0) - 500 Gpc-^r^^- The FOV of CTA 
with the shown sensitivity 5 deg would not be wide 
enough for a blind search. On the other hand, a survey 
mode with a wider FOV would not be sensitive enough to 
detect the signal. Thus, for CTA, a rapid follow-up ob- 
servation triggered by a wide-field X-ray telescope such 
as Swift or a Lobster-type instrument is needed. Assum- 
ing that the sky coverage is > 10%, one can expect > 0.2 
events yr~^ within 100 Mpc. The detection rate would be 
increased by a simultaneous operation of HAWC with a 



Fig. 3. — Energy fluxes of gamma rays corresponding to Fig[2] 
The emission duration is equal to that of X rays; t^y = 2.0 X 10^ sec. 
We show the cases of Dl = 10 and 100 Mpc (thick solid lines). For 
the former, we also show the injected spectrum without attenuation 
(dashed line) and only including attenuation within the emission 
region (thin line). The dotted lino shows 0.5 hr = 1.8 X 10'' sec 
differential sensitivity of CTA for a Scr detection. 

sensi tivity ^ 10~^°erg cm~^sec~^ for ~ 100 T eV gamma 
rays (|DeYoung fc HAWC CollaborationllMl) . 

A detection of the multi-TeV gamma-ray transient, as 
expected in this model, would also constrain the emis- 
sion mechanism of LL GRB. This is in contrast to the 
relativistic jet model, where as in the neutrino counter- 
part, the typical energy of the gamma rays injected by 
the photomeson reaction would be > PeV, which will 
be completely attenuated by the EEL even if they can 
escape the emission region. 

5. SUMMARY AND DISCUSSION 

We have shown that relativistic shock breakouts in SNe 
can be accompanied by multi-TeV neutrino and gamma- 
ray transients. These can provide diagnostics for a radia- 
tion mediated shock converting into a collisionlcss shock, 
and for baryon acceleration there. We can also get clues 
to the emission mechanism of LL GREs by detecting 
such high energy counterparts simultaneously with the 
prompt X-ray emission. The multi-TeV gamma rays can 
be detectable even from 100 Mpc away using CTA. These 
results motivate follow-up observations triggered by a 
wide-field X-ray telescope like Swift. 

While typically one expects very few neutrino events 
from those relativistic SNe, nevertheless searches for 
them would be aided by other possible counterparts. Us- 
ing the information of optical/infrared counterparts of 
core-collapse SNe, one can essentially fix the target po- 
sition within the angular resolution of IceCube/KMSnet 
< deg, and also restrict the time domain of the neutrino 
search within ~ day. The atmospheric neutrino back- 
ground (ANE) of IceCube/KMSnet within a circle of a 
degree over a day is roughly < 10^^ -Ei/,iooToV~^ events 
day~^. In terms of this ANE flux, neutrinos from rela- 
tivistic shock breakouts within Z^l '~ £^i/,iootcV Gpc can 
give a signal-to-noise ratio > 1 (see also Figl2]). One 
could then statistically extract 0(1) astrophysical neu- 
trinos by stacking the optical counterparts of O(IO^) 
SNe within z < 0.3, whether or not the X-ray coun- 
terparts are observed. Given that the whole sky event 
rate of such LL GREs would be 3 x 10^ yr"^ assum- 
ing Rhhiz = 0) ~ 500 Gpc^'^yr^^, a decadal SNe search 
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up to z < 0.3 with a sky coverage > 10% is needed 
While still a challenging task, this kind of a stronomy may 
be p ossible in the up-coming LSST era ()Lien fc Fieldsl 

Hooa). 

Ordinary SNe Ibc can also be accompanied by rela- 
tivistic shocks with /3sh ^ 0.1, which inevitably break 
out at a certain radius and produce neutrinos in a similar 
manner similar as discussed above. Although the typical 
fraction of the energy in the relativistic c omponent would 
be relatively smaU, £iso < 10^^ erg (e.g.. lSoderberg et al.l 



l2010f ). the event rate is much larger than for LL GRBs, 

Rihc{z = 0) ~ 2 X 10-* Gpc^Vr"^ (IMadau et al.l[T998h . 
Thus, they would give a comparable or a larger amounts 
of neutrinos by using the above stacking method. 

We thank Bin-Bin Zhang and Peter Veres for valuable 
discussions. This work is supported in part by a JSPS 
fellowship for research abroad, by NSF PHY-0757155 and 
the CCAPP workshop on "Revealing Deaths of Massive 
Stars with GeV-TeV Neutrinos" . 
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